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O ■ ABSTRACT 

We report on the identification of tlie liighest redshift submillimetre-selected source 
currently known: LESS J033229. 4— 275619. This source was detected in the Large Apex 
^ . Bolometer CAmera (LABOCA) Extended Chandra Deep Field South (ECDFS) Sub- 

5-H ' millimetre Survey (LESS), a sensitive 870-^m survey (crgro/jm ~ 1.2 mJy) of the full 

5^ ! 30' X 30' ECDFS with the LABOCA camera on the Atacama Pathfinder Experi- 

ment (APEX) telescope. The submillimetre emission is identified with a radio coun- 
terpart for which optical spectroscopy provides a redshift of z = 4.76. We show that 
the bolometric emission is dominated by a starburst with a star formation rate of 
^ lOOOAf0yr^^, although we also identify a moderate luminosity Active Galactic 
Nucleus (AGN) in this galaxy. Thus it has characteristics similar to those of z 2 
submillimetre galaxies (SMGs), with a mix of starburst and obscured AGN signatures. 
This demonstrates that ultraluminous starburst activity is not just restricted to the 
hosts of the most luminous (and hence rare) QSOs at z 5, but was also occurring 
in less extreme galaxies at a time when the Universe was less than 10 per cent of its 
current age. Assuming that we are seeing the major phase of star formation in this 
galaxy, then we demonstrate that it would be identified as a luminous distant red 
galaxy at z ^ 3 and that the current estimate of the space density of z > 4 SMGs 
is only sufficient to produce > 10 per cent of the luminous red galaxy population at 
these early times. However, this leaves open the possibility that some of these galaxies 
formed through less intense, but more extended star formation events. If the progeni- 
tors of all of the luminous red galaxies at z ~ 3 go through an ultraluminous starburst 
at z ^ 4 then the required volume density of z > 4 SMGs will exceed that predicted 
by current galaxy formation models by more than an order of magnitude. 

Key words: galaxies: high-redshift - galaxies: evolution - galaxies: formation - sub- 
millimetre - galaxies: individual LESS J033229.4-275619, CDS J033229.29-275619.5 
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1 INTRODUCTION 

It is now well established that a significant population of 
luminous red galaxies exist out to at least z ^ 1.5-2 and 
it has been claimed that most of these are massive and 
that some of them are old (> 10 ^ '^M(7) , ^ 2 Gyr) galaxies 



JCimatti et all 12004 [Paddi et all lioosl: iKong et alJ l2006l : 



iMcGrath. Stockton fc Canalizd l2007l : [Hartley et al.1 l20oi ) . 
There are indications that many of these massive galaxies 
were largely formed at earlier times, z ~ 2.5-3.5, with a frac- 
tion of them alread y having red colours, suggesting an aging 
stella r population ( Marchesini et al]l2007l : iMarchesini et al.l 
I2OO9I ). If correct, this suggests that these galax;ies must 
ha ve formed the bulk of their stellar populations at z > 
4 (iDaddi et al.l l2005l : ICimatti et al] l2008l : IStockton etHI 
I2OO8I ). The implied star formation rates (SFRs) for these 
progenitors would then be ~ 10^-10^ M0 yr~^ and if such 
vigorous activity occurred in a single halo, it is likely that 
these systems would appear as ultraluminous infrared galax- 
ies (ULIRGs), whose redshifted far-infrared emission would 
make them bright submillimetre galaxies (SMGs). However, 
the confirmation of large numbers of galaxies forming stars 
at a rate of ~ 10"^ M0 yr~^ at such high redshifts could pose 
a serious challenge to current hierarchical models, which pre- 
dict a z > 4 SMG surface density about an order of mag- 
nitude less (~ 10 deg~^) than would be implied if all mas- 
sive high-redshift galaxies form through this process (e.g . 



Granato et al . 2004; Baush ct al. 2005: H opkins et aLllioOsI : 
Bower et al ir2006: Swinbank et al. 2008 ) . 



Due to the negative K-corre ction, SMGs can e ssentially 
be detected out to 2 ~ 8 (see iBlain et al.ll2002l ). but the 
confirmation of SMGs (if they exist) at these redshifts poses 
a challenge. The large beamsize of submm cameras (e.g. full 
width at half maximum (FWHM) ~ 19" for LABOCA) 
precludes the determination of precise positions for large 
samples of SMGs without very signi ficant investment in 
sub- /millimetre interferometry (e.g. Downes et al.l 19991: 



Gear et al.ll2000l:lDannerbauer et al.i2002l . I2OO8I : IWang et al.1 



20071 : lYounger et al.l |2007| . l2008l )~This makes determining 
their nature through multiwavelength avenues extremely 
challenging, since there are typically many objects within 
this large beam that could be the source of the submillime- 
tre emission. The most successful route to identifying the 
counterparts of SMGs has been to search for 1.4-GHz coun- 
terparts (e.g. Ivison et al. 2002), exploiting the well-known 
far-in frared-radio correlation l|Condonlll992l : lHelou fc Bicavl 
1 19931 ). With deep, ^ 10 ^Jy rms, radio maps this provides 
precise positions for about 60 per cent of the SMG popula- 
tion {Ssso 4mJy), allowing follow-up spectroscopy and 
yielding a redshift distribution for radi o-identified SMGs 
peaking at 2 ~ 2.2 (jChapman et al.ll2005l ). With the advent 
of Spitzer, a similar approach has exploited 24 /xm counter- 
parts to SMGs, when radio counterparts were unavailable, 
incre asing the fraction o f iden t ified SMGs to ^ 70 per cent 
(e.g. Ilvison et al.1 12004 l2007l : iPope et al.1 12OO6I ). Unfortu- 
nately the radio and mid-infrared emission dim significantly 
with redshift, which means that the question of whether a 
significant proportion of SMGs lie at 2 ^ 4 remains open. 

Luminous submillimetre-emitt ing sources are of 
(e. 



course known at z 



4-6 



McMahon et al.l 



Omqnt et al.l Il996l: Ilvison et al.1 [l998bl: ICariUi et al 



2OOII : [Archibald et aLll200ll : iPriddev. Ivison fc IsaaklBool ). 



19941: 



200d. 



but these are all associated with rare and extreme Ac- 
tive Galactic Nuclei (AGN) whose space densities are so 
low that they are probably not a significant progenitor 
population for the luminous red galaxies seen at z ~ 3 
(although this conclusion depends upon the assumed duty 
cycle of the AGN). In contrast, the largest current redshift 
survey of th e more numerous subm illimetre-selected galaxy 
population jChapman et ahl |2005| ) identified no SMGs at 
z ^ 4. As noted above, the need for precise locations for 
the SMGs from their fijy radio emission biased the redshift 
distribution to z ^ 3.5 (|Chapman et al. 20051). and there 
could be a high-redsh i ft tail (iDannerbauer et al.l |2002|. 
20081 : IWang et al.l l2007l: lYounger et al.l l2007l: iGreye et al l 



20081: Iwaiig Barger fc Cowid l2009l: iKnudsen et al.l |2009| ) 



Pope et al.ri|2006h and lAretxaga et al.l l|2007h derived pho- 
tometric redshift estimates (albeit much less precise than 
spectroscopic redshifts) for possible counterparts to SMGs 
and proposed that a relatively low fraction (< 10 per cent) 
of SM Gs lie at 2: ^ 4, suggesting that the IChapman et al.l 
(|2005l ) redshift distribution is representative of all SMGs. 
Nevertheless, there has been recent con firmation of three 
radio-identified SMGs lying at z = 4-4.5 (jPaddi et al.i r2008l: 



ICapak et all I2OO8I : ISchinnerer et al.l I2OO8I ). indicating that 
a tail of SMGs does exist at z ^ 4. 

In this paper we focus on a luminous SMG discov- 
ered in our deep 870- /xm Large Apex BOlometer CAmera 
(LABOCA) Extended Chandra Deep Field South (ECDFS) 
Survey (LESS; Pis Small, Walter & Weiss; Weiss et al. in 
preparation, hereafter W09) which is associated with a faint 
radio and mid-infrared source. Spectroscopic and photomet- 
ric observations identify this source as an ultraluminous far- 
infrared galaxy ai z — 4.76. This is currently the highest 
redshift SMG known, suggesting that intense star formation 
was occurring in galaxies when the Universe was just over 
1 Gyr old. 

W e adopt cosm o logica l parameters from the WMAP 
fits in ISpergel etall (|2003D : = 0.73, = 0.27, and 
Ho = 71 kms~^ Mpc~^. At 2 = 4.76 the angular scale is 
6.5 kpc arcsec"^, and the Universe is 1.3 Gyr old. All quoted 
magnitudes are on the AB system unless otherwise noted. 



2 OBSERVATIONS AND SOURCE 

MULTIWAVELENGTH PROPERTIES 

The full area of the ECDFS was mapped at 870 fim with 
LABOC A on the 12-m A tacama Pathfinder Experiment 
(APEX: lGiisten et al.|[2006l ) telescope at Llano de Chajnan- 
tor in Chile as a joint ESO-Max Planck programme (the 
LESS collaboration). A total of 310 hrs of observing time 
(with 200 hrs of that on-source) was obtained which yielded 
a map covermg the central 30' x 30' of the field (0.25 deg^) 
to an average rms of ~ 1.2 mjy. This is the largest con- 
tiguous deep submillimetre survey yet undertaken (c.f. the 
~ 0.25 deg^ 2-mJy rms SCUBA Half D egree Extragalactic 
Survey - SHADES: [Coppin et al.ll2006l ). The description of 
the observations, reduction and analysis of the map is given 
in W09. The robust catalogue derived from this map con- 
tains 121 SMGs at > 3.7-(T significance (chosen to provide an 
estimated false detection rate of < 5 sources), and the mul- 
tiwavelength identification and characterisation of the SMG 
population in the ECDFS will be presented in subsequent 
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Figure 1. A selection of multiwavelength views of LESS J033229.4. The three panels show: (left) a true-colour ground-based BVR image 
of LESS J033229.4 from MUSYC, with a 2" x 2" inset illustrating the GOODS HST morphology in the F775W band (c onsistent with an 
unresolved point source); (centre) three-colour Spitzer 3.6, 4.5, 5.8 imaging from GOODS/SIMPLE with the lMiller e t al. ( 2008) VLA 
1.4-GHz radio contours overlaid (starting at 2 cr and increasing in steps of 1 it); (right) the 24 fira imaging from the Spitzer-FYDFih survey 
with the LABOCA 870/im emission from W09 overlaid as contours (starting at 2 a and increasing in steps of Icr). The submillimetre 
position is indicated by a cross, and we also plot a circle representing the 8" search radius used to locate potential counterparts. We have 
robustly identified the SMG with a radio and 24-^m counterpart ~ 1.5" away, which is coincident with a. z = 4.76 galaxy (see text). 
Each panel is 15" X 15" in size and is orientated with North to the top and East to the left. 



papers. Here we focus on the multiwavelength characteris- 
tics of a single source: LESS J033229.4-275619 (hereafter 
LESS J033229.4), detected at 870 with a flux density of 
6.3±1.2mJy (S/N = 5.1) at 03''32"29.41= -27°56'18.90" 
(J2000), with a 90 per cent positional uncertainty of 6" 
(WOO). WOO apply a flux-deboosting algorithm to correct 
the raw flux and noise values for all of the LESS sources in 
ord er to obtain the mo st accurate submm photometry (see 
e.g. lCoppin et al.ll2005l ). In this paper we use the deboosted 
flux and noise of 5.0 ± 1.4 mjy for LESS J033229.4. 

We use a new radio catalogue (Biggs et al. in prepara- 
tion) which comprises emitters with peak flux densities in 
excess of 3 a, where a is determined locally, ex tracted from 
the V ery Large Array (VLA) 1.4-GHz map of iMiller et all 
l|2008h . We search within 8" of the submm position for 
potential radio co unterparts (see e.g. iDownes et al.lll986l : 
llvison et al.ll200"3 ) and identify a single candidate at a ra- 
dial distance of 1.5" with a flux density of 24.0 ± 6.3 /xJy 
(S/N=3.8), assuming that it is a point source. The proba- 
bility of this radio emitter being aligned by chance with the 
SMG is less than 2 per cent and we can assign this radio 
counterpart to the SMG with some confidence. Radio im- 
ages suffer from similar flux boosting effects as the submm 
images and we thus adopt a deboosted 1.4-GHz flux density 
of 18.8 ± 6.3 /^Jy from Biggs et al. (in preparatio n) , follow- 
ing ex tensive simulations of the kind described bv llbar et al.l 
l|2009l) . The positional accuracy of the radio data allows 
us to pinpoint LESS J033229.4 to qj2ooo = 03''32™29.30" 
5j20oo = -27°56'l9.40 " with an uncerta inty of Aa = 0.2" 
and AS = 0.3" (see llvison et all |2007| and Fig. [ij. Us- 
ing this precise position for the robust radio counterpart 
of LESS J033229.4 we can now use the impressive archival 
observations of the ECDFS to derive its spectral energy dis- 
tribution. 

In the mid-to-far-infrared, a counterpart is de- 
tected at 3. 6-8Atm using the InfraRed Array Cam- 
era (IRAC; iFazio et al] |2004| ) data from the Great 



Observat ories Origins Deep Surve y (GOOD 



Dickinson. Giavalisco fc the GOODS teami \200. . 

Giavalisco et al.l |2004| ) imaging of the central part of 



the ECDFS (GOODS-South) and in the wider area cover- 
age provided by the Spitzer IRAC/MUSYC Pubhc Legacy 
Survey in the ECDFS (SIMPLeQ; Damen et al. in prepara- 
tion) . For the present analysis we use the 3-" radius aperture 
photometry in the SIMPLE catalogue, as it provides close 
to optimal S/N for point sources, with the appropriate 
aperture corrections applied and corrected to total fluxes. 
In the mid-infrared, LESS J033229.4 is detected by the 
Far-Infrared Deep Extragalactic Legacy survey (FIDEL; 
Dickinson et al. in preparation) i n the Multiband Imaging 
Photometer for Spitzer (MIPS; iRieke et al.l l2004h 24 fim 
catalogue at 03''32™29.30' -27°56' 18.78" , with a positional 
accuracy of 1-2", given the FWHM and the source's S/N. 
The MIPS source lies within 0.6" of the radio counterpart 
and hence the radio and MIPS emission originate from the 
same source, given the combined error circle. We report a 
24^im flux density of 824^,^ = 31.6 ± 5.1/iJy (S/N = 6.6), 
exploiting the deeper GOODS imaging (Chary et al. 

in preparation). LESS J033229.4 is not detected at either 
70 or 160 fim in the deep FIDEL imaging and so we adopt 
3-0- upper limits at these wavelengths. 

At shorter wavelengths, the region around 
LESS J033229.4 is covered by optical and near-infrared 
imaging from several ground and space-based public sur- 
veys: the Mult i wavel ength Survey by Yale-Chile (MUSYCH; 
iGawiser et al] I2OO6I ). the Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) Galaxy Evolution 
from Morphologies and SEDs survey (GEM^l iRix et all 



http: / / www. stsci . edu / science /goods/ 

http: / /data. spitzer. caltech.edu/popular/SIMPLE 

http:/ /www. astro. yale.edu/MUSYC 

http: / / www.mpia.de / GEMS / gems. htm 
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bOO^ I. and GOODS (|Giavalisco et alj |2004| : Giavalisco et 
al. in preparation). We list the fluxes or limits from these 
surveys in Table [1] The source is undetected in all bands 
shortward of ~ 600 nm, although it is significantly detected 
in the redder wavebands, Veoe, ^^,1775, zgso. Unfortunately 
the source falls ~ 1" outside the deep Very Large Telescope 
(VLT) Infrared Spectrometer And Array Camera (ISAAC) 
JHK imaging of GOODS-South and so there are currently 
only wea k constraints on its near-infrared fiuxes from 
MUSYC (|Tavlor et al.ll2008l ). 

Finally, turning to the X-ray waveban d: from their 
analysis of the 1-Ms Chandra image, Giacconi et al.1 (j2002l l 
report the detection of an X-ray source at 03''32™29.44' 
-27°56"20.18" with a positional uncertainty of 2.1", with 
0.5-2 and 2-10 keV fluxes of (1.6 ± 0.5) x 10"^® and < 7.6 x 
10~^^ ergs~^ cm~'^, respectively. This position is consistent 
with the GOODS optical counterpart 0.8" away. However, 
this source (ID 618) is only listed in th eir supplementary 
catal ogue and is based on a SExtractor l|Bertin fc Arnoritg 
Il996l ) analysis of the X-ray image, for which the software was 
not intended. In contrast, no X-ray source is detected at this 
position in any of the analyses of the 250 ks , 1-, and 2- Ms im- 
ages o f th e CDFS, bylLehnier et al.1 (|200i),|AlexandereEaL| 
l|2003h or ILuo et alj ( 20081 '). respectivelv. We thus extract 
3- ff upper limits f rom the deepest 2-Ms Chandra imaging 
of ILuo et ahl (|2008l ). corresponding to rest-frame 3-11.5 keV 
and 11.5 — 46keV luminosity limits of < 2.1 x 10"'^ and 
< 1.5 X lO^^ergs"^, respectively, which we use to constrain 
any contribution from an AGN to the bolometric luminosity 
(see ^ [3:211 . 

Table [T] summarises the fluxes or limits derived from 
the multiwavelength coverage of LESS J033229.4. From the 
optical photometry it is clear that the source is a V- 
band dropout. Assuming this arises from the presence of 
the Lyman limit in the IZ-band, this suggests a redshift 
of z ~ 5. The morphology of the source in the redder 
HST imaging is very compact, with an observed FWHM 
of % 0.1" (or a half-light radius of ^ 0.3 kpc), indicat- 
ing that the source is unresolved (see inset in Fig. [T|. In- 
deed, the source was photometrically and morphologically 
selecte d as a candidate z = 3.5-5.2 AGN bv lFontanot et al] 
(|2007l ). On this basis, LESS J033229.4 was observed with the 
FOcal Reducer and low dispersion Spectro graph (FQRS2; 
lAppenzeller et al.l Il998h on the VL T by IVanzella et"al] 
l|2006t ) (see also IVanzefla et al.l l2008h as part of the ESO 
Large Programme 170.A-0788 (PI: Cesarsky). Their ID for 
this source is GDS J033229. 29-275619.5 and they obtained 
a 4-hr exposure of it using the 3001 grism, deriving a redshift 
oi z = 4.76 on the basis of the identification of strong Lya 
emission, Nv and a continuum break (Fig. [2]). The photo- 
metric properties of LESS J033229.4 are consistent with this 
claimed redshift, in particular the discontinuity at ~ 600 nm 
is well reproduced by the Lyman limit at z — 4.76 (Fig. [3]). 
This makes this source the highest redshift submillimetre- 
selected SMG yet discovered and indeed one of the highest 
redshift 24-/im-detected sources so far identified. 

We obtained additional spectroscopic observations of 
LESS J033229.4 using the DEep Iinaging Muhi-Object Spec- 
trograph (DEIMOS; iFaber et al] 120031 ') on Keck on the 
nights of 2007 October 10 and 11. The total exposure time 
was 6.5 hrs in good conditions using the 600 lines mm~^ 
grating and the resulting spectrum (Fig. [2]) shows contin- 



Table 1. The observed photometry for LESS J033229.4 from Chan- 
dra, HST, MUSYC, Spitzer, APEX and the VLA. Note that we list 
the deboosted submm and radio fiuxes, which is important in order 
to correctly model the SED and investigate the radio and submm 
propertie s of the SMG (see text). We use the HST- ACS photome- 
try from IStark et al.l | |2007|) , which has already been corrected for 
the Lya emission measured from the FORS2 rest-frame UV spec- 
trum. We also give the 3.6 /im flux density, corrected down by 30 per 
cent to reflect the typical equivalent wid th of H« in z ~ 2 SMGs 
jSwinbank et al.ll2004l : iTakata et al.ll2006l) . Where the source is un- 
detected we quote a 3-0" upper limit to the flux density. Note: 1/jJy 
is the equivalent of 23.9 AB magnitudes. 

Filter Flux Density Reference 



/Wavelength 


(M Jy) 


0.5-8 keV 


< 3.7 X lO"^** erj 


0.5-2 keV 


< 8.5 X 10^1^ erj 


2-8 keV 


< 6.3 X 10"^** erj 


U (352 nm) 


< 0.087 


S435 (433 nm) 


< 0.013 


B (461 nm) 


< 0.039 


V (538 nm) 


< 0.059 


Veoe (597 nm) 


0.072 ± 0.006 


R (652 nm) 


0.215 ± 0.020 


»775 (771 nm) 


0.350 ± 0.013 


Z850 (905 nm) 


0.405 ± 0.015 


J (1.2/jm) 


< 1.65 


if (2.lMm) 


< 3.20 


3.6 fim 


2.86 ± 0.06 


4.5 ^m 


4.04 ± 0.08 


5.8 /im 


6.3 ± 0.4 


8.0 fim 


9.2 ± 0.4 


24 ^m 


32 ± 5 


70 fjm 


< 2500 


160 /im 


< 33000 


870 Mm 


5000 ± 1400 


20 cm 


18.8 ± 6.3 



' cm^'' Luo et al. (20081 
' cm^^ Luo et al. (2008) 
' cm~^ Luo et al. (2008) 

Gawiser et al. (_2006) 

^bntano^^t_al^ (2007) 
Gawiser et al. (2006) 

^amise^_e^_al^ 1^2006) 

^tar^^^^ ^2007) 

^awise^_e^_al^ ||2006) 

^tar^^^^ (2007) 
Star k et al. (2007) 
Tayl or et al. (2008) 
Taylor et al. (20081 
Damen et al. in preparation 
Damen et al. in preparation 
Damen et al. in preparation 
Damen et al. in preparation 
Chary et al. in preparation 
Dickinson et al. in preparation 
Dickinson et al. in preparation 
W09 

Biggs et al. in preparation 
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Figure 2. The optical spectrum of LESS J033229.4 derived from 
combining our new dee p DEIMOS sp e ctrum with the archival 
FORS2 spectrum from iVanzella et al.l | |2006|) . The upper spec- 
trum is the combined dataset at the FORS2 resolution (~ 
1.3 nm), while the lower spectra (arbitrarily offset in flux for 
clarity) show the individual DEIMOS spectrum convolved to 
~ 0.8 nm resolution and the FORS2 spectrum. The combined 
spectrum displays two strong emission lines and an associated 
continuum break. We identify the lines as Lya and Nv, yielding 
a redshift of 2 ~ 4.76. We mark these two emission features on 
the spectra and also indicate the expected wavelengths for other 
potential emission and absorption lines at the adopted redshift 
oi z = 4.76: Siiil26.0, Oi 130.3, Ciil33.6, Siiv 140.3, Feni 149.4, 
Sill 152.7 and Civ 155.1, as well as atmospheric absorption at 
760 nm (shaded bar) . The upper scale gives the rest-frame wave- 
length assuming a redshift of z = 4.76, and the flux scale is nom- 
inally in ergcm"'^ s~^ A~^. 
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uum emission and the two strong emission lines seen in the 
FORS2 spectrum. We r etriev ed the reduced FORS2 spec- 
trum of IVanzeUa et alJ (j2006l ) from the ESO archive and 
then convolved the DEIMOS spectrum to the ~ 1.3 nm reso- 
lution of the FORS2 spectrum before averaging the two spec- 
tra. The combined spectrum has an effective 10-m equiv- 
alent integration time of ~ 9hrs and is shown in Fig. (2] 
The combined spectrum displays a narrow Lya emission 
line at 700.5 nm and a broad Nv emission line at 714.5 nm 
with a width of ~ 5nm or ~ 2000 km s~^ in the rest- 
frame and an integrated flux density comparable to Lya 
(~ 1.4 X 10~^^ ergs~^ cm~^). To search for weak broad lines 
we have convolved the spectrum to the FWHM of the Nv 
line and find an emission feature corresponding to Civ 155.0 
at very low significance, likely due to the strong telluric 
emission at these wavelengths. The spectral properties of 
the source, in particular the relatively strong and moder- 
ately broad Nv, indicates the presence of an AGN, though 
the weak X-ray and 24 fim flux densities relative to the op- 
tical and far-infrared bands suggest a relatively weak AGN 
contribution to the bolometric luminosity. 

High-redshift AGN ar e often associated with rest-frame 
far-infrared emission fe.g. Ilvison et al]|2008l ). and as a final 
test of our identification of the SMG with this source, we can 
also ask what is the likelihood of finding such a high-redshift 
source close to the submillimetre position by chance . Usin g 
the AGN luminosity function from Fontanot et al.l l|2007h . 
we calculate how many 2 > 4 AGN counterparts we would 
have expected to find at random given their volume den- 
sity and the offset of the submillimetre source to the AGN, 
e = 1.5". LESS J033229.4 has a rest-frame 145-nm lu r ninos - 
ity of M145 = -21.5, calculated following [Fan et all (|2000l ') 
and adopting fi, oc i^"" with ct = 0.5. The volume den- 
sity of AGN brighter than LESS J033229.4 is ^(Mws < 
—21.5) ~ 2.6 X 10"®Mpc~'^, equivalent to a surface den- 
sity N ~ 34deg~^ for the integrated density at z > 4. The 
Poisson random probability of finding a z > 4 AGN at least 
as bright as ours is P = 1 — exp[— vr^^A^] ~ 10~^. Therefore 
it is extremely unlikely that our SMG is located within 1.5" 
of a high-redshift source by random chance, confirming the 
proposed association of the submillimetre emission with the 
z = 4.76 source. 



3 RESULTS & ANALYSIS 

3.1 Spectral Energy Distribution 

The bright submillimetre fiux density of LESS J033229.4 
suggests that a massive burst of star formation is under- 
way in this system ai z — 4.76. To test if the observed radio 
and far-infrared emission is consistent with being produced 
by star formation, we calculate the radio-to-submm spectral 
index for LESS J033229.4, deriving afi' = 1.01 ± 0.11. This 
spectral index is in good agreement with the mean spec- 
tral index a t z = 4.76 predicted from 17 local star-forming 
galaxies by ICarifli fc YunI l|2000l ): af^f = 1.04 ± 0.13, as 
weU as with local ULIRGs such as Arp220, aff ~ 1.09. 
This confirms that the far-infrared-radio emission from 
LESS J033229.4 is consistent with massive star formation. 

We therefore plot the multiwavelength spectral energy 
distribution (SED) of LESS J033229.4 in Fig. |3] and com- 
pare this to the SEDs of Arp220 and M82, normalised 



to match the observed SsTOMm of LESS J033229.4 - corre- 
sponding to rest-frame 150-/im emission, near the peak of 
the dust emission. This shows that while the Arp 220 and 
M82 SEDs both fit the radio-far- infrared emission from 
LESS J033229.4, M82 is significantly brighter in the mid- 
infrared, with only Arp 220 providing an adequate descrip- 
tion of the SED across the rest-frame optical-to-radio. This 
confirms that this z — 4.76 SMG has an SED consistent with 
local star-formation-dominated ULIRGs, although we note 
that the rest-frame UV and 24 /im emission exceed that seen 
for Arp 220. Therefore we have also overplotted a redshifted 
composite QSO SE P from the atlas of local, optically bright 
quasars compiled bv lElvis et al.l (|l994l ). scaled to match our 
observed rest-frame UV photometry, in order to assess the 
likely bolometric contribution of an unobscured AGN to our 
observed SED (see Fig. [3]). It is apparent from this sim- 
ple exercise that any contribution from an unobscured AGN 
to the infrared luminosity of LESS J033229.4 is likely to be 
negligible, and we investigate this in detail in ii l3.2l 

From the scaled Arp 220 template, we calculate Lfir — 
^{8-1000 Mm) = 6.1 X 10^^ Lp,, a factor ~ 5x the bolomet 



ric lumino sity of Arp 220 (jSoifer et al.l Il984l : I Sanders et al.l 
Il988l . l2003l '). We note that adopting the scaled M 82 template 
results in an Lpir a factor of ~ 1.5 x higher. Owing to the 
negative K correction in the submillimetre waveband, the 
bolometric luminosity of LESS J033229.4 is clos e to the me- 
dian l uminosity o f the z ^ 2 populat i on in .Chapman et ahl 
(|2005l ) (see also iKovacs et all I2OO6I : ICoppin et al.l |2008| ). 
even though its luminosity distance is 2.5 x larger. 

Based on these SED arguments, which indicate that 
much of the far-infrared emission arises from star formation, 
we derive a SFR of ^ lO OOMnyr"^ for LESS J033229.4, 
following iKenniciittl lll998h who ass umes a burst lifetime 
of ~ 100 Myr and a ISalpetej (jl955l ) initial mass function 
(IMF). Next, we estimate a dust mass of Md ~ 5 x 10** Mq 
using the 870 fim fiux density and the usual relation S'i/„t.. = 
B^, (T )k^, Md (1 + z)/Df^ (e.g. lHuehes. Dunlop fc Rawlingd 
Il997f ). where Dl is the cosmological luminosity distance 
and B^i is the Planck function evaluated at the emit- 
ted frequency, v' = fohs (1 -I- z). We have extrapolated 
from an average wavelength-dependent mass-absorption 
coeff icient of ki25„^ = 2.64m^k g~^ assuming jS — 
1.5 (jDunne. Eales fc Edmuridsl l2003l ). Based on the sim- 
ilar far-infrared luminosity and spectral properties of 
LESS J033229.4 to z - 2 SMGs we adopt the average 
gas/dust ratio for this population (~ 60: ICreve et a 1120051 : 
iTacconi et aLll2006l : iKovacs et allbooel : ICoppin et alll2008l ). 



and so predict a gas mass for LESS J033229.4 of ~ 3 x 
10^" M0 . Comparing this mass to the SFR, we expect that 
LESS J033229.4 would consume its gas at its present SFR 
in ~ 30-40 Myr. This predicted gas mass is similar to 
that determined for other z > 4 SMGs ()Daddi et al.ll2008l : 
[Schinnoror ct al. 2001). 

To place limits on the potential stellar mass of the 
galaxy, we now investigate if the rest-frame UV-optical 
emission can be fit by a stellar population usin g the Hy- 
perZ package /Bolzone lla. Miralles fc Pefldl2000l ). We stress 
that this is a necessary, but not a sufficient, requirement 
to show that this emission actually arises from stars, rather 
than the AGN. We assume a solar metallic ity stellar popula- 
tion model from iBruzual fc CharlotI (|l993l 'l. a lCalzetti et al.l 
(|2000l ) starburst attenuation law, and either a constant or 



6 Coppin et al. 




Observed X [ |im] 

Figure 3. The SED of LESS J033229.4 showing the dominant and luminous rest-frame far-infrared emissi on indicating that the source 
is an ultraluminous infrared galaxy. For comparison we have overplotted the SED of Arp220 and M82 from lSilva et al l Il998^ red s hifted 
to 2: = 4.76 and scaled to match the Sgro^m flux density of LESS J033229.4 and also a composite QSO SED from lElvis et al.l ill 9941 ') 
redshifted to z = 4.76 scaled to match the rest-frame UV flux density. The good agreement with the starburst templates in the far- 
infrared and radio suggest that the bulk of the bolometric emission in this system arises from star formation. We also plot the SED 
of the z = 4.5 SMG SMM J100054.5-I-023436 from Capak et al. (2008; COS) to show the similar far-infrared properties, but the much 
more intense rest-frame UV emission from that system. The inset (in the same units as the larger plot) shows the best-fitting stellar 
SED from HyperZ to the optical— 8 /^m photometry, a 40-Myr old constant star formation model with moderate reddening. Ay ~ 1.5, 
which provides an acceptable fit of the rest-frame UV-optical emission in LESS J033229.4. However, we caution that the current SED 
can be equally well-fit by a power-law, as expected from a reddened AGN, leaving open the possibility that a significant fraction of the 
rest- frame near-infrared emission arises from the AGN. 



single-burst star formation history. Under these assump- 
tions, HyperZ shows that the observed SED can be fit 
by a stellar population a.t z = 4.76, although statistically 
a reddened power-law provides as good a fit. The best-fit 
SED is a 40-Myr old constant star formation model with 
moderate reddening, Ay ~ 1.5 (see inset Fig. [Sjl. Ifowever, 
the absence of useful limits on the galaxy's SED in the ob- 
served JHK bands, sampling the critical region around the 
Balmer and 4000A breaks, means that there is significant 
uncertainty in these parameters, with statistically accept- 
able solutions with ages up to ^ 500 Myr and lower red- 
dening, Av ^ 0.5. Even this range of fits should be used 
with caution as the potentially dominant contribution from 
the AGN, as well as the complex mix of dust and stars, in 
this SMG will mean that the luminosity-weighted ages and 
reddening will vary strongly with wavelength as different re- 
gions within the galaxy are probed. Hence, in the following 
we simply ass ume a canonical dura tion of the SMG-phase 
of ~ 100 Myr (|Swinbank et al.ll2006l ). and focus on the rest- 
frame near-infrared luminosity which should suffer the least 
AGN contamination and provide the most robust estimate 
of the stellar mass of the galaxy. 

From a simple interpolation of the SED we derive a 
rest-frame K-hand luminosity of Mk ~ —24.0 or Lk ~ 
5 X 10^^ Lq. To estimate a limit on the stellar mass from 
this we use the parameterisation for Lj^/Mp for a con- 
stant star formation model from iBorvs et al.l (|2005l ). For 



star formation which has been occuring for 100 Myr this 
gives Lk / M* 10 (40- and 300-Myr durations correspond 
to La'/Mq ~ 50 and 7), implying a typical stellar mass of 
M* 5 X 10^" Mq . This has an uncertainty of at least a fac- 
tor of 3-5 X, even before considering the potential contribu- 
tion from an AGN to the rest-frame near-infrared emission. 
Nevertheless, we note that this limit on the stellar mass is 
comparable to the predicted gas mass of the system and that 
for ages of 100 Myr, the stellar mass produced in the burst 
implies a SFR of ~ 500 Mq yr ~^, similar to that d erived from 
the far-infrared. In contrast, IStark et al.l (|2007| ) estimate a 
stellar mass for LESS J033229.4 of M* ~ 1.3 x lO" Mq, as- 
suming the mass-to- light ratio appropriate for a ~ 1-Gyr old 
stellar population, which is unlikely to be valid for a system 
with a current SFR of - 1000 Mq yr"^ 

3.2 Constraints on AGN energetics 

The rest- frame UV compact morphology (FWHM of <, 0.1"; 
or a half-light radius of ^ 0.3 kpc) and spectroscopy of 
LESS J033229.4 indicate that it hosts an AGN. For com- 
parison, about ~ 10 per cent of z ~ 2 SMGs arc cla ssed as 
'com pact', with half-light radii of < 0.6 kpc (Swinban k et al.l 
l2009h . Similar spectral signatures of AGN emission are 
also seen in the rest-frame UV spectra of around 25 per 
cent of the radio-det ected SMG population at 2 2-3 
(iChapman et aLll2005h . suggesting a mix of star formation 
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and AGN in the galaxy. In particular, LESS J033229.4's 
spectral properties are very similar to the first spectro- 
scopically identified SMG: the type-2 or narrow-line QSO 
SMMJ02399-0136 at 2 = 2.803 jlvison et al.1 Il998al ). 
SMM J02399— 0136 also shows strong Lya and a comparably 
bright and broad Nv line, with a width of 1800 km s~^. Al- 
though the emission lines in SMM J02399— 0136 are propor- 
tionally stronger relative to the continuum and the source 
is intrinsically more luminous than LESS J033229.4. 

We can use the strength of the Nv emission and other 
SED constraints to estimate the potential luminosity of the 
AGN component in LESS J033229.4 and compare these to 
our X-ray limits, where we are confident that any AGN 
emission will dominate. Using the observed 2-8 keV limit 
(rest-frame 11.5-46 keV) and assuming F = 2.0 we derive 
a conservative limit on the rest-frame 2-10 keV luminosity 
of < 1.7 X 10*'*ergs~^. This will be almost unaffected by 
the presence of absorption assuming absorbing column den- 
sities of A'h % 10^^ cm~^, i.e. not heavily Compton thick. 
A more sensitive limit comes from the observed 0.5-2 keV 
limit which yields a rest-frame 2-10 keV luminosity limit of 
< 2.5 X lO^'^ergs"^. However, this is also more sensitive to 
the assumed absorption in the system and so we will use our 
more conservative limit in the following. 

Starting from the measured Nv emission-line luminos- 
ity of ~ 3 X lO*'^ ergs~^, we can use the observed Nv-X-ray 
luminosity ratio found fo r SMM J02399-0136 (.Ivison et all 
ll998al : lBautz et al.llioool ') to predict the rest-ft-ame 2-10 keV 
luminosity of LESS J033229.4 of ~ 1 x 10*" ergs'^ We can 
also attempt to estimate the luminosity of the AGN using 
the 24 fim flux density under the extreme assumption that 
the 24 fim emission is entirely due to the AGN. Using the de- 
rived rest-frame 6 fj,m luminosity of 1.2 x 10*^ ergs~^ we es- 
timate a rest-frame 2-10 ke V luminosity of « 3 x 10*" erg s~^ 
based on [Alexander et al. 

I |2008). Similarly, assuming that 
the UV continuum luminosity of LESS J033229.4 is dom- 
inated by emission from the AGN, we can predict the 
AGN luminosity. We take the observed rest-frame 145- 
nm absolute magnitude (M145 — —21.5) and convert this 
to rest-frame 250-nm a ssuming a = —0.44 {Si, cx z/"; 
IVanden Berk et al.lbOOll '). We then use the Qox relationship 
of lSteffen et al.l "( 20061 1 and convert between different X-ray 
bands assuming F = 2.0, to derive a rest-frame 2-10 keV 
luminosity of 5 x 10*^ ergs~^. 

These three estimates are in the range 0.5-3 x 
10** erg s"^ and so are consistent with our observed X-ray 
limit, given the significant uncertainties in the calculations. 



4 DISCUSSION 

Our estimates of the stellar and gas masses of 
LESS J033229.4 in §0 are both crude and highly un- 
certain; nevertheless, both masses are comparable to those 
derived for the better-constrained and similarly bolomet- 
rically luminous 2 ~ 2 SMG population. The implied 
baryonic mass of the galaxy is thus likely to be ^ 10^^ Mq, 
suggesting that the ultraluminous starburst is occuring in 
a massive galaxy. Indeed, similar constraints on the host 
galaxy masses ha ve been derived for the other exarnples o f 
2 > 4 SMGs (e.g. lDaddi eraLll2008l : ISchinnerer et al.ll2008h . 
Our three predictions of the X-ray luminosity of the 



AGN in §[33] suggest the luminosity of the AGN in 
LESS J033229.4 is ~ 10**ergs-^ and that the AGN may 
dominate the emission in the X-ray and UV, while making 
a modest contribution to the mid-infrared luminosity. Using 
the X-ray-to-infrared luminosity ratio for local AGN (e.g. 
[Alexander et al.ll2005l ) we expect that the AGN contributes 
^ 20 per cent of the bolometric luminosity and hence the 
bulk of the bolometric luminosity of LESS J033229.4 arises 
from star formation. 

The discovery of LESS J033229 . 4 at z = 4. 76 and three 
z = 4-4.5 SMGs from'Capa k et al.1 (|2008l ) and lDaddi et al.1 
(2008) indicates that galaxies with SFRs of 1000 Mq yr"^ 
are present at high redshifts. Using these four z > 4 SMGs, 
and the areas of the surveys they were derived from, we 
can place a lower limit on the surface density of SMGs 
at these redshifts of ^ 6deg~^, or a volume density of 
^ 1.5 X 10~^ Mpc"'' conservatively assuming the surveys are 
uniformly sensitive across 2 = 4-8. Note that if we include 
the likely 2 > 4 SMG froin IWang et af] (|2007| ) (see also 
IWang. Barger fc Cowidl2009l ). then these limits increase by 
~ 15 per cent. 

To test whether these high-redshift SMGs can evolve 
into luminous red galaxies at 2 ~ 3, we compute the ex- 
pected fading of LESS J033229.4 from 2 ~ 4.5 to 2 ~ 3. 
Qualitatively, the inferred SFRs of the four 2 > 4 SMGs 
(~ 10^ M0 yr~^), with estimated duration of order 100 Myr, 
are sufficient to form > 10^^ Mq of stellar mass needed to 
match the properties of luminous red galaxies seen at 2 ~ 3, 
and there is enough time between 2 ~ 4 and 2 ~ 3 to age 
their stellar populations sufficiently. To quantitatively con- 
firm this we adopt an approach that, while uncertain, is 
insensitive to the AGN contribution to the rest-frame K- 
band luminosity of LESS J033229.4. We assume a constant 
star formation rate at the observed level of 1000 Mq yr~^ for 
100 Myr and that the starburst is observed half way through, 
then we estimate fading after ~ 1-1.5 Gyr of AK ~ 1.1- 
1.3 in the rest-frame if-band resulting in a typical SMG- 
descendent at 2 ~ 3 with A'Ik ~ —24.7 ±0.1 and observed 
colours of (J — K) ~ 0.6-0.9, assuming Bruzua l fc CharlotI 
(1993) stellar models. Using the .Marasto n ( 1 9981 ) stellar 
population models these predictions become AK ~ 0.7- 
1.7, Mk -24.7 ± 0.5 and {J - K) ~ 1.0-1.3. These pre- 
dicted luminosities and colours are in reasonable agreement 
with those of the reddest luminous galaxies at 2 ~ 2.5-3.5, 
which have Mk J=, - 23.5 ± 0.5 and {J - K) > 1.3 (e.g. 
iMarchesini et al.ll2009h . 

We can also ask whether the space densities of these two 
populations are consistent with this evolutionary cycle. Our 
lower limit on the space density of SMGs at 2 ~ 4-8 is ^ 
1.5 X 10~^ Mpc""^, so correcting for the duty cycle assuming 
a 100- Myr lifetime (~ 9x) this suggests that there will be 
^ 10~^ Mpc~^ descendents at 2 ~ 2.5-3.5. This compares 
to the measured space density of luminous galaxies at these 
redshifts of ~ 5 x lO"'^ Mpc~^ of which ~ 20-30 per cent 
have red colours (|Marchesini et al ] l2007l : IMarchesini et al.l 
l2009h . or ~ lO"'^ Mpc"^. This suggests that we have so far 
only detected ^ 10 per cent of the SMG progenitors at 2 > 4 
of the 2^3 evolved luminous galaxies if all of the latte r 
have a previous ultraluminous phase (c.f. IStark et al.ll2009l ). 
or alternatively that there should be a surface density of 
~ lOOdeg"^ SMGs at 2 > 4, corresponding to ~ 10 per cent 
of the ^ 5mJy SMG population. 
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Figure 4. A comparison of the volume density of radio-detected 
SMGs with 5i.4GHz > 30 ^Jy and Sss Oum > 5mjy from the 
GALFORM model of lBaueh et alj ||2005| ') to our hmit on the vol- 
ume densit y of 2: > 4 SM Gs and the measured density of 2 ~ 2 
SMGs from lChapman et al. (2005). This demonstrates that these 
detailed theoretical models are just consistent with these current 
observations, and so the discovery of more examples of 2 > 4 
SMGs within the ECDFS would result in a disagreement with 
the model predictions. 



More detailed theoretical galaxy formation models have 
achieved a rough agreement with the number count s and the 
median redshift of the SM G po pulation (e.g. Baugh et al.l 
I2OO5I : ISwinbank et al.|[20oi ). The lBaugh et al.l (|2005l ) semi- 
analytical galaxy formation model implements a top-heavy 
IMF in bursts in order to reprodu ce the submillimetre nu m- 
ber counts and has been shown bv lSwinbank et ah I dlOOS') to 
provide a reasonable description of the far-infrared proper- 
ties of submillimetre galaxies at z ~ 2. We now investigate if 
the proportion of spectroscopically-confirmed high-redshift 
SMGs can provide a more stringent test of this model. In 
Fig.|4]we compare the volume densities of spec troscopically- 
confi r med SMGs to the model predictions (|Baugh et al.l 
l2005l '). lBaugh et all (|2005l ) predict surface and volume densi- 
ties of ^ 7deg"^ and ^ 2x 10"^ Mpc"^ (for z_ ~ 4-8), respec- 



tively , for SMGs with Sgso nm > 5 mjy (see ISwinbank et al.l 
l2008f ). Imposing a radio-detection constraint (^ 30 /^Jy), 
which can be directly compared against the SMG samples 
with spectroscopic confirmation, reduces these predictions 
by a factor of ~ 2 (Fig.|4]). Both of these estimates are con- 
sistent with our current limits on the volume densities of 
SMGs at z > 4 (as well as the me asured volume density at 
z ^ 2 from IChapman eral]l2005l ). However, we note that 
the foUow-up of the LABOCA survey of the ECDFS has 
only just begun and so it is possible that more 2: > 4 SMGs 
will be uncovered. Hence, we suggest that it is probable that 
the current semi-analytic models underpredict the numbers 
of very high redshift ULIRGs. In particular, if an example 
of a z > 5 SMG is found in existing submillimetre surveys 
then this will place strong constraints on these models. 



5 CONCLUSIONS 

LESS J033229.4 is currently the highest redshift 
submillimetre-selected galaxy known. At z = 4.76, it is seen 
just over 1 Gyr after the Big Bang, yet it has similar proper- 
ties to the SMG population seen at the peak of their activity 
at z ~ 2 when the Universe was nearly 3x older. The SED 
and rest-frame UV spectrum indicate that LESS J033229.4 
is a composite AGN/starburst galaxy caught in a phase 
where star formation (SFR~ lOOOMoyr"^) is dominating 
the bolometric emission. 

At z — 4.76, 870 /im samples near the SED peak, pro- 
viding a good constraint on the far-infrared luminosity of 
SMG of 6 X lO^'^ L0 . But in order to accurately deter- 
mine the SED peak shape and bulk dust temperature of 
LESS J033229.4 one requires photometry on the Wien side 
of the peak of the SED from observations of the ECDFS 
at 250, 350 or 500 /im with the Herschel Space Observatory 
or the Balloon -borne Large Apert ure Submillimetre Tele- 
scope (BLAST: Fascale et al.ll2008l ). Our preferred Arp220- 
like SED fit predicts 250, 350 and 500 /im flux densities of 
~ 5, 9, and lOmJy at these wavelengths, which are likely to 
prove to be challenging for either Herschel or BLAST due 
to confusion. Although a deep 350 /im (rest-frame 60/im) or 
450/im flux density measurement would be feasible with fu- 
ture APEX measurements using the SABOCA submillime- 
tre camera or the James Clerk Maxwell Telescope's Sub- 
millimetre Common-User Bolometer Array 2 (SCUBA-2), 
respectively. 

Based on our estimated dust mass of 5 X 10* Mq and 
adopting a 2: ~ 2 SMG average gas/dust ratio, we pre- 
dict a gas mass of 3 x 10^° Mp), similar t o that measured 
for o ther z > 4 SMGs l|Daddi et al.ll2008l : ISchinnerer et all 
l2008h . Atz = 4.76 ^^CO(5-4) emission should be detectable 
using current interferometric facilities. The confirmation of 
a massive reservoir of gas in LESS J033229.4 would enable 
the study of the gas and dynamics of galaxy formation in 
detail when the Universe was only 1 Gyr old, and so better 
place LESS J033229.4 in context with SMGs and red lumi- 
nous galaxies at z ~ 3. 

As with the original identification of SMGs 
jlvison et al.l Il998al ). it is likely that selection effects 
may be strongly influencing our view of the properties 
of the very high-redshift SMG population. We note that 
if LESS J033229.4 had not been identified as a compact, 
potential high-redshift, AGN then it would not have been 
spectroscopically observed and if it had not shown strong 
line emission (at least in part from the AGN) then a redshift 
would never have been measured. Simil arly, the morphologi - 
cally complex optical counterpart of the lCapak et al.l (|2008l ) 
SMG was first identified as a weak radio source associated 
with a potentially z > 4 y-ba nd dropout o b ject, and the 
redshifts for the two z = 4.05 iDaddi et af] |2008l ) SMGs 
were obtained from the serendipitous detection of redshifted 
CO (4-3) emission in observations of an unrelated nearby 
z — 1.5 BzK galaxy. A complete view of the properties 
of the z > 4 SMG population will require a concerted 
identification programme and multiwavelength follow-up, 
but not withstanding this, at present we can conclude 
that there appears to be an similarly diverse population of 
intense star-bursting galaxies at z ^ 4 as there is at z ~ 2. 

Assuming that red luminous galaxies at 2 3 are all 
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descendants of z = 4-8 SMGs we estimate the space den- 
sity of the SMG progenitor population would be of the order 
of ~ 10~^ Mpc~^ at 2; ~ 4-8, assuming a 100-Myr lifetime 
for the ULIRG-phase, corresponding to a surface density of 
~ 50-100 deg"'^. Current observational limits on the z > A 
SMG population are thus consistent with them being the 
progenitors of > 10% of the passive massive galaxy pop- 
ulation at 2 ~ 3. These limits are in agreement with the 
semi-analytical model predictions for radio-detected SMGs 
at these redshifts. However, we note that if further follow-up 
shows that the surface density of SMGs at z ~ 4-8 is closer 
to the ~ 100 deg~^ necessary for them to form all of the red 
luminous galaxy population at z ~ 3, then the models will 
be underpredicting the number of SMGs at z > 4 by about 
a order of magnitude. Alternatively, many z ~ 3 red lu- 
minous galaxies may form without a superluminous phase. 
Either way, statistically significant samples of SMGs need 
to be studied and followed up with spectroscopy in fields 
with extensive multiwavelength coverage before we will be 
able to determine true predominance and role of SMGs at 
z > 4. The panoramic and deep Cosmology Legacy Surveys 
planned with SCUBA-2 will obtain the tens of thousands of 
SMGs needed to begin the search in earnest for SMGs at 
even higher redshifts, z 3> 4-8. 
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